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Abstract. We consider discrete-time nearest-neighbor quantum walks on random environments 
in one dimension. Using the method based on a path counting, we present both quenched and 
annealed weak limit theorems for the quantum walk. 1 1 1 1 

1 Introduction 

Classical random walks have proved to be a very important tool for various fields. Quan- 
tum walks are the quantum counterpart of classical random walks. There are two types of 
quantum walks. One is the discrete-time walk and the other is the continuous-time one. Re- 
views and books on quantum walks are Kempe [I] , Kendon [2] , Venegas-Andraca [3] , Konno 
in [5] , for examples. Ambainis et al. [6] intensively studied discrete-time quantum walks on 
Z, where Z is the set of integers. Here we consider a discrete-time nearest-neighbor quan- 
tum walk in a random environment (QWRE) on Z. We present weak limit theorems for the 
QWRE, both quenched (i.e., conditional upon the environment) and annealed (i.e., averaged 
over the environment), by counting the number of paths that takes a quantum walker from 
the origin to a position. Numerical results on continuous-time QWRE were reported in Yin 
et al. [7]. As for the classical case, the study of random walks in random environments was 
initiated in the mid-1970s and has been investigated by many researchers. Sinai's random 
walk in a random environment [S] is one of the typical models. See Zeitouni [H] for a review 
of random walks on random environments. To the best of our knowledge, any rigorous result 
for the QWRE is not known. Therefore we hope that our investigation would be a first step 
towards the study of the QWRE. 

The rest of the paper is organized as follows. Section 2 gives the definition of the QWRE. 
In Sect. 3, we present the quenched weak limit theorem (Theorem 13. 2p . As a consequence, 
the annealed weak limit theorem (Corollary 13. 3p is also obtained. Section 4 is devoted to 
the proof of Theorem 13.21 

Abbr. title: Quantum walks on random environments 

AMS 2000 subject classifications: 60F05, 60G50, 82B41, 81Q99 

PACS: 03.67.Lx, 05.40.Fb, 02.50.Cw 

Keywords: Quantum walk, random environment, limit theorem, Hadamard walk 



1 



2 Definition 



The time evolution of the QWRE on Z is determined by a sequence of 2 x 2 random unitary 
matrices, {Ux x E Z,} , where 



^x dx 



with ax,bx,Cx,dx € C and C is the set of complex numbers. The subscript x indicates the 
location. The unitarity of Ux gives 

\(^x\ ~l~ \^x\ \bx\ ~l~ Mx'l 1; ^x^x ~\~ bxdx 0, Cx ^x^x) dx ^x^xy 

where z is the complex conjugate of z E C and A^. = det Ux = cbxdx — bxCx with | A^-l = 1. 

The discrete-time quantum walk is a quantum version of the classical random walk with 
additional degree of freedom called chirality. We should remark that there is a strong struc- 
tural similarity between quantum walks and correlated random walks, see Konno [10]. The 
chirality takes values left and right, and it means the direction of the motion of the particle. 
At each time step, if the particle has the left chirality, it moves one step to the left, and if 
it has the right chirality, it moves one step to the right. Let define 



1^) 
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\R) 
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so Ur acts on two chiralities as follows: 



Ux\L) = ax\L) + Cx\R), Ux\R) = bx\L) + dx\R), 

where L and R refer to the left and right chirality state, respectively. To define the dynamics 
of the QWRE, we divide Ux into two matrices: 



ax bx 




Qx 





C.7: dx 



with Ux = Px + Qx- The important point is that Px (resp. Qx) represents that the particle 
moves to the left (resp. right) at position x at each time step. 

We will explain the QWRE more precisely. Let fl = M^, where R is the set of real 
numbers. A random environment is an fi-valued random variable u = {Ux : x E 1^} with 
probability measure P. We will assume that P is a product measure on (fi, J-"), where J-" 
is the Borel cr-field of Q. So we write P = Ylxez^^' where Px is a probability measure on 
(R, Q) and Q is the Borel cr-field of R. Throughout this paper, we consider the following Ux'- 

A typical example is that u = {ux : x G Z} is independent and identically distributed. That 
is, Px does not depend on position x. In particular, when a; = 0, i.e., Ux = {x E Z) or 
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Px = Sq, where da is the delta measure at a e M, the dynamics is non-random and position- 
independent. This case is equivalent to the Hadamard walk determined by the Hadamard 
gate Ux = H: 

1 1 1 
1 -1 



V2 



The walk has been extensively investigated in the study of the quantum walk. 
The set of initial qubit states at the origin for the QWRE is given by 

$={v. = ^[«,/3]eC^:|ap+|/3p = l}, 

where T is the transposed operator. In the present paper, we take = ^[l/\/2, i/\/2] as 
the initial qubit state of the QWRE. Then the probability distribution of the Hadamard 
walk starting from at the origin is symmetric. We call the quantum walk the symmetric 
Hadamard walk here. 

Let En{l,iTL) denote the sum of all paths starting from the origin in the trajectory con- 
sisting of / steps left and m steps right. In fact, for time n = l + m and position x = —l-\-m, 
we have 

^nKh in) — ^x(h)^x{mi)^x(l2)^x{m2) " " " -^a;(i„_i)^a;(m„_i)-^0 Vo ' 

where the summation is taken over all Ij, rrij > satisfying /i -|- ■■■ + /« = /, mi -|- • • • -|- m„ = 
m, Ij + rrij = 1, and \x{aj+i) — x{bj)\ = 1, x{aj+i), x{bj) EZ{a,b E {l,m}). We should note 
that the definition gives 

H„+i(/, m) = Px+l Enil - 1, m) + Qa:-1 En{l, m - 1). 

For example, in the case of / = 3, m = 1, wc have 

54(3, 1) = Q^sP^2P-iPo + P-1Q-2P-1P0 + P-iPoQ-lPo + P-iPoPiQo- (2.1) 

The probability that a quantum walker in an environment G is in position x at time 
n starting from the origin with (p*{— '^[l/^/2,i/^/2]) e $ is defined by 

Pni^n = x) := ||S„(Z,m)(/7*||^ 

where n — l-\-m and x — —I -\- m. 

For each a; e fi, P!^ is a probability measure on (Z, where "H is the set of all subsets 
of Z. Thus we can define a probability measure P„ :— P ® P^ on O x Z by the formula 

¥^{H X {X, e A}) := / P^{{X^ e A}) P{duj), 

for any H E H, A G 1j. Statements involving P^ are called quenched. On the other hand, 
statements involving P„ are called annealed. Expectations under P and P^ are denoted by 
E and E!^, respectively. In n = 4 case, by definition, a direct computation implies 

Pr(x. = -4) = i±9M, jy(x, = -2) = ^±i|HM, 
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If the probability measure of Uq, i.e., Pq, is symmetric, (for example, ojq is a uniform distri- 
bution on [— 7r,7r]), then 

P4(X4 = -4) = 1, P4(X4 = -2) = ^, P4(X4 = 0) = |, 
P4(X4 = 2) = ^, P4(X4 = 4) = 1. 

In fact. Corollary 14.41 tells us that Pn(X„ = ■) is the same as the probability distribution of 
the symmetric Hadamard walk for any n > 0. 



3 Results 

In this section, we will present both quenched and annealed weak limit theorems for the 
QWRE. Quantum walks behave quite differently from classical random walks. For example, 
in the classical case, the probability distribution is a binomial one. On the other hand, the 
probability distribution of the Hadamard walk has a complicated and oscillatory form. For 
the classical case, the well-known central limit theorem holds. That is, y/n converges 
to a normal distribution as n — > oo, where Yn denotes the probability distribution of the 
location for the random walker at time n. Concerning the quantum walk defined by 2 x 2 
non-random unitary matrix Ux = U, Konno [HI |12j showed the corresponding weak limit 
theorem. In particular, for the symmetric Hadamard walk, i.e., a; = 0, we have 



Theorem 3.1 



lim P°(u< — <v]= [ fxix) dx, 

n^oo \ n J 



where — oo < u < v < oo, and 

Here Ia{x) = 1 (x e A), = {x ^ A). 

We can extend this result {u = 0) to the general u & Q. The following quenched weak limit 
theorem is our main result of this paper. 

Theorem 3.2 (Quenched case for Xn/n). For any fixed environment u &Vt, 
lim \u < — - < v] = [ {1 — sin(a;o)a;} fxix) dx, 

n^oo \ n J 



where — oo < u < v < oo. 

One of the interesting points of this result is that the limit density depends only on uq. By 
Theorem 13. 2[ we can obtain an annealed weak limit theorem immediately as follows: 
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Corollary 3.3 (Annealed case for Xn/n) 



lim P„ ( n < — < ) = / {1 - E{sm{uJo))x} dx, 



n 



where — oo < u < v < oo. In particular, when the probability measure of uq, (i.e., Pq) is 
symmetric, we have 



lim P„ ( u < — - < v] = [ fxix) dx. 
n^oo \ n J 



If the probability measure of uq is symmetric, then the limit density is the same as that of 
Theorem 13.11 for the symmetric Hadamard walk. 
If Ux = Ti/2 {x E Z), then 

In this non-random case, when we take = '^[1/ \/2,i/\/2] as the initial qubit state, the 
weak limit theorem shown by Konno [HI [T2] gives the same conclusion as Theorem 13. 2[ 
That is, the limit density is {1 — x) fxix). 

From Theorem 13. 2[ for each environment w G f2, the variance of the limit density, V^, is 
given by 



Corollary 3.4 



2-V2 / 2-V2 . 
V = < 1 sm (woj 



When Uq = (i.e., Pq = 5q), takes the maximum value. Therefore the additional noise 
makes the variance smaller. On the other hand, if ujq = 7r/2 (i.e., Pq = then 
takes the minimum value. So the additional noise makes the variance larger. As for random 
walks in random environments, for example, the following result for the Sinai random walk 
corresponding to a usual symmetric random walk (i.e., in a non-random environment) is 
known. Let denote the probability distribution of the location for the random walker 
at time in a random environment. Then Z„/(logn)^ converges to some non-degenerate 
random variable as n — )■ oo, (see [H |9] for more details). The result tells us that this order 
(logra)^ is smaller than the order ^/n for the non-random environment case. On the other 
hand, our Theorems 13.11 and 13 . 2 1 imply that both quantum walks in non-random and random 
environments have the same order n. The difference between them is the magnitude of the 
variance as shown in Corollary 13.41 Moreover, compared with the non-random environment 
case, the variance for the corresponding random environment case can be larger or smaller. 

In our setting, {Ux ■ x E is a. sequence of position-dependent 2x2 random unitary 
matrices. On the other hand, when {f/„ : n = 1, 2, . . .} is a sequence of time- dependent 2x2 
random unitary matrices, Mackay et al. [13] and Ribeiro et al. p3] investigated and showed 
that the probability distribution of the quantum walk converges to a binomial distribution 
by averaging over many trials by numerical simulations. Konno [15] proved their results by 
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using a path counting method. In this paper, we use a similar approach to prove Theorem 

We discuss other related works [HI [13 UHl [IS [201 EI] • Brun et al. [HI HI] considered 
discrete-time quantum walks with many coins on the line. By a different method, they 
analytically showed that the number of different coins has one of the same effects noted in 
our work: the more different coins are used, the smaller the variance, while still retaining 
the quantum scaling. The corresponding weak limit theorem was presented in [TH]. Richter 
[T^ [20] proved that decoherent quantum walks on finite sets of positions under repeated 
randomized measurements have the potential speed up a large class of their classical coun- 
terparts. This work is still analytic using yet another method. Ermann et al. [21] investigated 
discrete-time quantum walks on the line in chaotic and regular environments by using the 
method given in [Tni[IT]. The family of chaotic environments includes the so-called quantum 
multibaker map as a special case. Their analytical results suggest that chaotic environments 
induce decoherence more effectively than regular ones. 



4 Proof of Theorem 3.2 



In this section we will give the proof of Theorem 13.21 First we introduce the useful random 
matrices to compute S„(/,m): 







Sx 





^x bx 



In general, we obtain the following table of products of the matrices Px,Qx,Rx and 
(x G Z): for any x, y G Z, 





Py 


Qy 


Ry 


Sy 


Px 


OixPy 


^xRy 


Ry 


bxPy 


Qx 


CxSy 


dxQy 


^xQy 


^xSy 


Rx 


CxPy 


dxRy 


^x Ry 


^xRy 


Sx 


OixSy 


^xQy 


^xQy 


bxSy 



where PxQy = bxRy, for example. From this table and Eq. (12. ip . we obtain 

24(3, 1) = (6_ic_2a-i + a-i&oC-i)Po + ct-iao&i-Ro + c_3a_2a-i5'o. (4.2) 

We should note that Po,Qo,Ro and Sq form an orthonormal basis of the vector space of 
complex 2x2 matrices with respect to the trace inner product {A\B) = tT{A*B), where * 
means the adjoint operator. So S„(/,m) has the following form: 

E:„(/, m) = pn{l, m)Po + qn{l, m)Qo + r„(Z, m)Ro + Sn{l, m)So. 

As in a similar way for the one-dimensional discrete-time quantum walk (see Konno, [ITl [T2] , 
for example), we can compute explicit forms of p„(Z, m), g„(/, m), r„(Z, m) and s„(/,m) as 
follows: 
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Proposition 4.1 When I Am(:= min{/,m}) > 1, we obtain 

p^(l,m)^e'^\l,m)pl^\l,m), qr,(l,m) ^ el^\l,m) qi''\l,m), 
rn(l,m) = eW(Z,m) rf )(Z,m), Sn(l,m) = el:\l,m) sl''\l,m), 



where 



and 



oii^-i+'-'+'^-ii-m-i)) ifl^i>rn, 

e-i(^o+-+^rn-i) ifl-l<m 
e'('^°+-+^-(i-m)) ifl>m-l, 



/ \ n-l (/-l)Am , . , . 



7=1 
-1 iA(m-l) 



7-1/ 

^i'"('-)K7l)"""^S"<-^'-'C"-i)("o")- 

When I Am — 0, we have 

pi^)(0,n) = rf H0,n) = siPiO,n) = 0, )(0,n) = (-l/v^)"-\ 
e(f)(0,n) = e([)(0,n) = e^^)(0,n) = 1, e(,5)(0,n) = e-^('^i+-+'^"-i), 

anc? 

p(f )(n, 0) = {l/V2r-\ £\n, 0) = rW(n, 0) = )(n, 0) = 0, 
e(f)(n,0) = e-(--+-+--("-)), ei^)(n,0) = eM(n,0) = e(:)(n,0) = 1. 

We should remark that if a; = (i.e., Hadamard walk), then Q^^(Z,m) = Qn\l,'rn) — 
Qn\l,rn) — Qn\l,m) = 1. Therefore we use superscript (H) for a^\l,m) (a = p,q,r,s). 
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For example, in / = 3,m = 1 case, Proposition 14.11 gives 

ef (3, 1) = e'"-S (3, 1) = e^("°+^-i+"-2), 

Q^;\3, 1) = e^(^°+"'-^\ el'^(3, 1) = e*('^-i+'^-2), 

p^(3,l) = 2(l/y2)^ gf)(3,l) = 0, rf ^(3, 1) = ^(3, 1) = 

On the other hand, by Eg. (14.21) . 

24(3, 1) = {1/V2f (2e''"-iFo + e'("°+"-i)i?o + e'^'^^'^''-'^ So) . 
So we have the same conclusion. From Proposition 14. we have 
Corollary 4.2 Ifn = l + m and x = —I + m, then 

P„-(X. = x) = P°(X„ = x) + i - H/,^)'} ■ sin(a;o) 

ri^^(/,m) 



+q[^\l,m) 53 (^(1 + e^^'o) B^^'^H^H ei''^(/,m; 
=: P0(X„ = x) + Wr{uJo) + 1^2(0;), 
where '^{z) is the imaginary part of z EC 

From now on, we will show that W2{oj) = for any l,m. First we consider I — 1 > m. In 
this case, we have 

e W (/, m) = e^('^o+-+'^-('-™-i)) , e^f ) (/, m) = e^('^-i+-+'^-('-™)) . 
By using the above equations. 



53 



Jl + e^^"") e(f)(/,m) ei''^(/,m)j = 53 ((1 + e^*"") Qi'\l,m) el?\l,m) 

= $5 ((1 + e2*^o)e-^'^°) = 0. 

Therefore we obtain W2{uj) = 0. Next we consider I — 1 = m case. Then we have 

e(f)(/,m) = 1, Ql'\l,m) = e^('^o+^-i), 
eM(/,m) = e^-«, eW(Z,m) = e^--. 

These imply 



53 



(1 + e^^"") e(f)(^,"^) Qn\l,m)) = 53 ((1 + e^*"") G^'^H^,^) 0n^(^m) 



= $5 ((1 + e2*'^o)e-^'^°) = 0. 

So we get W2{uj) = 0. For other cases, we have the same conclusion in a similar fashion. As 
a consequence, 
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Proposition 4.3 

P^iX^ = x) = P°(X„ = x) + l {pi!'\l,mr - )(/,m)2} ■ sinM, 

where n = I + m and x = —I + m. 
Then we have immediately 

Corollary 4.4 If the probability measure of Uq (i.e., Pq), is symmetric, then 

P„(X„ = x) = P°(X„ = x), 

for any n > and x (zl^. 

Let [x\ denote the integer part of a; G M. For 1 <l <[n/2\, pi^\l,m) and q^\l,m) can 
be rewritten as 

n-l I 



( H ] ( H ] 

Furthermore we will rewrite pn {l,m) and {l,m) by using Jacobi polynomials. Let 
P^'^{x) denote the Jacobi polynomial which is orthogonal on [—1, 1] with respect to (1 — 
xY{l + x)^ with jy, n > —1. Then the following relation holds: 

P^'^ix) = r(!i+|)r("^ + i) n + u + ^^ + l;u+l;{l- x)/2), (4.3) 

where 2Fi{a, b; c; z) is the hypergeometric series and T{z) is the gamma function. In general, 
as for orthogonal polynomials, see [22]. Then we see that 

= il/2)-^'-'\F,i~{k - 1), n - A; + 1; 2; 1/2) 
= ^(l/2)-('=-i)p«(0). (4.4) 



7=1 
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The first equality is given by the definition of the hypergeometric series. The second equality 
comes from the following relation: 

2Fi(a,6; c; z) = {1 - 2;)-%Fi(a, c - 6; c; z/{z - 1)). 

The last equality follows from Eq. (I4.3p . In a similar way, we have 

E (-1)^-^ I C 7 - ^ = (l/2)-(-^)p°--"(0). (4.5) 

By using Eqs. fl4.4p and (14.51) . we get 
Lemma 4.5 For 1 <l < [n/2], 



^) (-ir'{p.°:r'(o)-P/_:r'(o)} 



^ \ n-2«+l 



-ir" {(V)i^.^T"(o)-P.°:r"(o)}. 



Finally we will consider the limit behavior of the characteristic function E^(e^^^"^'^) with 
^ G M as n — i- oo. From Proposition 14. 3[ we have 

E-(e^«^"/'^) = E°(e*«^"/") + e^C{n-2/)/n (pf _ if _ gf - /)2} . 

1=0 

(4.6) 

For the first term of the right-hand side of Eq. (14. 6 p corresponding to the symmetric 
Hadamard walk, Theorem 13.11 implies 

lim E°(e*«^"/") = / e*«" dx. (4.7) 

"^"^ J-i/y/2 7r(l - a;2)A/l - 2x2 

Concerning the second term of the right-hand side of Eq. (14. 6p . by a similar argument in 
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[ni[T2], we see that 



1 " 
lim - > 

n-5>oo 2 ^ 
1=0 

n 

= lim y 



Mn^2i)/n [piH){i^n- if - qf\l,n- if] 



1=0 



(l+l/v^)/2 



(l-l/V2)/2 



X 



(l+l/V2)/2 

(l-l/V2)/2 
1/V2 



J^{n-2l)/n 



n — I 



n-2l+2 r 



pl,n.-2i 
M-1 



r)pl,n-2l p0,n-2l . ( p0.n-2l 



,iai-2x) 



pl,n-2l 
M-1 



2 ^ 2{n — I) i^n-2l p0,n~2l . ( p0,n-2l^ ^ 



2/ TTy/l - 2(1 - 2x)^ 



1 — X 



X 



J({l-2x)_ 



2x 

1 — X 

1 



2 1 



1 — X 



X 7 2(1 -x) 



dx 



2x - 1 



3 _ 



-1/V2 7r(l-x2)Vl-2x2 



27r^l - 2(1 - 2x)2 x(l 
dx, 



(ix 



X 



X 



(4. 



where P/l" = P/l" ^'(0) {i = 0, 1). The first equahty comes from Lemma [4.51 Combining 
Eqs. dMD and (giD with Eq. (gS]) gives 

1/V2 



lim E;^(e^«^"/") : 

n— ^oo 

SO the proof of Theorem 13.21 is complete. 



■.^x 1 - sin(u;o)x 
-i/v^^ 7r(l-x2)Vl-2x2 ' 
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